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Abstract: 

We study the clustering of colour strings in the transverse plane of nucleus-nucleus 
£ — | collisions and argue that deconfinement sets in when the string density within a cluster 

reaches a critical value. We show that this implies a minimal cluster size at the onset of 
^ ■ deconfinement, which in turn results in discontinuous behaviour for J/i/j suppression. 
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Statistical quantum chromodynamics predicts that strongly interacting matter will 
become deconfined at high temperatures and/or densities. The aim of high energy nuclear 



p , 

collisions is to produce and study this new state of matter, the quark-gluon plasma, 
in the laboratory. One prerequisite for such studies is an understanding of the onset 
of deconfinement. In the most general sense, the deconfinement of quarks and gluons 
means that these partons are no longer constrained to the distributions by which they 
are governed within hadrons [p]]. Such constraints are presumably removed once enough 
nucleon-nucleon interactions overlap in space and time, i.e., for a sufficiently high collision 
density. A partitioning of the coloured partons into colour-neutral subsets is then no 
longer meaningful; instead, there emerge colour-conducting clusters of much larger than 
hadronic size. 

Cluster formation is the central topic of percolation theory, and hence deconfinement 
and percolation appear to be closely related 0-0]. The aim of the present note is to 
pursue this relation in more detail, in particular for the finite-size environment encountered 
in nucleus-nucleus interactions, and to apply the results to recent data on J/ ip suppression 
in Pb — Pb interactions || . 

Percolation theory || in its conventional continuum form [|7| investigates the formation 
of clusters of differents sizes, if N extended objects ('discs' or 'spheres' in two or three 
space dimensions, respectively) are randomly distributed in a given spatial region, allowing 
overlap. For reasons to be explained shortly, we consider here the two-dimensional case, 
taking the spatial region to be a circle A = irR 2 of radius R and the discs as circles 
a = irr 2 of radii r (see Fig. 1). The basic question is how the average geometric size A c \ of 
a cluster of connected discs varies with the overall disc density N/A\ it is convenient here 
to use the dimensionless density n = a(N/A). The percolation point n p is defined to be 
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that density n at which in the limit R 



n < n p ) or last time (when n > n p ); in two dimensions, it is found to be n p 



oo infinite clusters appear for the first (when 

1.175 §, 

in accord with numerical results [J7|. Approaching n p from below, the average cluster size 
thus diverges, 
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while coming from above, the ratio of average cluster size to overall area vanishes, with 
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n > n p ; 



(2) 



T\ and r 2 denote the relevant critical exponents. In all previous work on percolation 
in strong interaction thermodynamics it was assumed that deconfinement sets in 

at the percolation point of some specific extended objects. In general, this conjecture 
does not seem tenable; in the Ising model, for example, random percolation and thermal 
critical behaviour (spontaneous magnetization) do not coincide || ||, and pure SU(N) 

The reason for 
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gauge theory is in the same universality class as Z N spin systems 
the difference lies in the different concept of cluster in the two cases, which tends to 
force thermal critical behaviour to higher densities than required for percolation. Hence a 
redefinition of percolation clusters is needed to make the two phenomena coincide; on the 
lattice, suitably defined site-bond clusters are shown to do this [[|. The generalization 
of such considerations to the percolation of extended objects is not clear, however, so 
that deconfinement could presumably set in at n p or at some higher density; for further 
Here we therefore leave the value of the deconfinement threshold n r 



discussion, see 
open. 
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Figure 1: Clustering of discs in a circular 



Consider now a central high energy collision of two identical heavy nuclei, during 
which essentially all nucleons undergo several interactions. These occur in such rapid 



succession that the nucleons cannot 'recover' between successive interactions. Each in- 
dividual nucleon-nucleon collision establishes a colour flux tube or string between the 



collision partners; the additive quark model (T^j suggests that at present (SPS) energies 



it connects two triplet colour charges. The transverse radius r of such strings is expected 
to have the form [Vj, lj 

r 2 = — ln(L/L ), (3) 
no 

where L is the separation between the colour charges and L < 0.1 — 0.3 fm the 'string 
formation' length JT5|, IB ; a ~ 0.16 GeV 2 denotes the string tension. Extensive lattice 



QCD studies find that in the range 0.5 < L < 2 fm, r ~ 0.2 — 0.3 fm and varies at 
most weakly pj |TB[. In a dense environment, colour screening will moreover constrain 



the possible values of L, so that here a radius of the mentioned size appears appropriate. 
For nuclear collisions in their early stage, we thus obtain a spaghetti-like structure of 
intertwined QCD strings, and a cut in the transverse plane results in the picture shown 
in Fig. 1, giving a distribution of transverse string areas over a circle of nuclear area. We 
want to study what happens when these overlap more and more to form connected spatial 
regions of increasing size in the transverse planeQ; hence the appropriate framework is the 
two-dimensional case of percolation theory introduced above. 

We shall now first study how the cluster density in the transverse plane varies as 
function of the density of discs, then how the geometric cluster size is related to the 
cluster density. It is useful to define the dimensionless cluster density 

n cl =(^)a> 1, (4) 
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where N c \ denotes the average number of discs making up the cluster of size A c \, by 
definition, n c \ = 1 for an isolated disc. Similarly, we define dimensionless cluster size 
measures 

s c i = — > 1 (5) 
a 

and 

Sci = ^ < i; (6) 

these are normalized such that s c i = 1 for an isolated disc and S c \ = 1 when the cluster 
covers the whole basis area A. To simplify notation, we shall from now on always mean 
the average over many configurations when we speak of quantities such as density, cluster 
density or cluster size. 

We shall see shortly that the cluster density n c \ increases smoothly with the overall 
density n. Since also the cluster size s c i grows with n, clusters of a certain density will 
always have a certain size. Hence at the onset of deconfinement, at some n c \ = n% clusters 
will have an intrinsic size, s c cl = s c i(n^). Deconfinement therefore begins in regions which 
cannot be arbitrarily small; the new phase makes its first appearance in an area of finite 
size. This is the main result of our note and leads to interesting consequences for the 
pattern of J/ip suppression by deconfinement. We now turn to cluster formation in detail. 

To consider a specific case, we set R/r = 20, corresponding to R = 5 fm and r = 0.25 
fm. The area of the discs thus is 1/400 of the 'nuclear' cross section A = nR 2 over which 

lr The first study of this kind was carried out for a specific phenomenological string model [0. 



they are distributed; we shall return later to the role of R/r and its specific value. To 
determine average cluster quantities, we calculate the size of the cluster containing the 
origin of A (the shaded area in Fig. 1) in a given configuration and then average over many 
configurations ||. In Fig. 2a we see that the cluster density n c \ increases monotonically 
with the overall density n, starting from unity for n — > 0. 

The crucial result of percolation theory is that the size s c \ of the cluster, in contrast to 
its density, shows a dramatic variation with n. In particular, for R/r — ► oo, s c \ diverges 
at the percolation point n p = 1.175, i.e., when the overall density is somewhat above one 
disc per disc area; S c \ vanishes there. Of course there is a considerable overlap of discs, 
and at n p , the ratio of the area covered by discs (in clusters of all sizes) to the overall area 
becomes 1 — exp[— 1.12] ~ 0.67; it converges to unity only for n — > oo. For finite R/r, this 
critical behaviour will be modified by finite size effects, and for the value R/r = 20 chosen 
above, we obtain the cluster size behaviour shown in Fig. 2b. As expected, it shows the 
strongest variation around the percolation point; for sufficiently large systems it is here 
governed by Eq. (|l|) for n below and Eq. (||) above n p . 
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Figure 2: The cluster density n c \ (a) and the cluster size Sci('Sci) (b), as functions of the 
overall density n, for R=5 fm, r=0.25 fm and a uniform random distribution. 
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Let us now try to 'translate' n and n c \ into variables relevant for deconfinement in 
nuclear collisions. The number of strings is determined by the number of nucleon-nucleon 
collisions; this can be determined experimentally by measuring the rate of Drell-Yan dilep- 
ton production. Nucleon-nucleon collisions and strings play a basic role in the early stages 
of the nucleus-nucleus interaction; the modification of the parton distributions associated 
to the onset of deconfinement can occur through successive nucleon-nucleon interactions, 
if the time between the collisions is too short for the restoration of a physical nucleon 
Eventually the strings fuse into 'wounded' nucleons which lead to the observed soft 



hadron production, and on this level the number of times a nucleon was hit by others 



no longer matters. In a Glauber-based description [IS], both the average (transverse) 
collision density and the corresponding density of wounded nucleons can be calculated 
for nucleus-nucleus collisions at fixed impact parameter b. Through convolution with the 
correlation between the number of wounded nucleons at fixed b and the energy of 
the measured secondary hadrons, it then becomes possible to determine all quantities of 
interest in terms of experimental observables. Fig. 2a thus effectively shows the collision 
density in the average cluster as function of the overall collision density; the latter can 
then be correlated to the density of wounded nucleons and thus to E T . 

As already noted, deconfinement is expected to set in when there is enough 'internet- 



ting' between interacting nucleons [T8|], i.e., when the collision density within a cluster 
becomes sufficiently high; we denote this value by n c d . It is reached at a certain value 
n c = n(n c cl ) of the overall density, which can, but does not need to be the percolation 
point (see Fig. 2a). By Fig. 2b, this identifies the corresponding critical cluster size 
s ci = s ci(^d)- m other words, requiring a critical collision density automatically forces 
the cluster to have a certain minimum size. This result clarifies the question of how much 
medium is needed before one can speak about a new phase: systems of geometric size 
s c \ < s c cl will on the average not have a string density sufficient for deconfinement. 

The onset of deconfinement in nuclear collisions is thus governed by two parameters: 
the deconfinement string density n c cl (or some equivalent quantity) and the transverse 
string size r; with n c cl and r fixed, we can then determine the size s c cX of the bubbles 
of deconfined medium present at the deconfinement threshold. The size of the nuclear 
transverse area irR 2 is given by nuclear geometry. For a realistic description, we use 
Wood-Saxon nuclear distributions [H7H and a Glauber-based formalism to calculate the 



distribution of collisions in the transverse plane as function of the impact parameter in 
A — B collisions [T?| . The distribution of the string cross sections in the transverse plane 
is thus no longer completely random; it is weighted with the mentioned nucleon-nucleon 
collision density. 

In the remainder of this note, we shall apply our considerations to the study of J / ip sup- 
pression [19| in nuclear collisions as observed in S — U and Pb — Pb data 0. For this 
purpose, we will choose r in the range 0.2 - 0.3 fm and fix n c cl such that central S — U 
collisions are just below the deconfinement threshold [2IJ. We can then address two 
questions: 

• How is the deconfinement threshold thus defined related to the percolation point? 

• What is the effect of the finite cluster size on the deconfinement pattern? 

As answer to the first question, in previous work |4]] the coincidence of deconfinement 
and percolation had been assumed, n(n^) = n p . On the second question, a recent study 
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20| had already shown that a minimal bubble size for the deconfmed medium leads to 
an abrupt onset of J/ip suppression as function of the collision centrality. The critical 
bubble size was there obtained through superheating in a first order phase transition. 
Such a picture is, however, not as easily adapted to the environment produced in nuclear 
collisions as the present cluster formation approach, which moreover does not contain the 
bubble size as an open parameter. 



Our application will be based on the successive suppression scenario [£2|] - f24| , in which 
the x an d the ip' are dissociated essentially at the deconfinement point, the J/tp later. 
This implies that the fraction (about 40 %) of J/^'s coming from x and ip' decay will be 
suppressed when n c \ = n% if these states are produced within a deconfmed cluster. The 
directly produced J/^'s then survive until n^ lss ( J/ij))] based on energy density estimates 
p2| , p3f] , this is expected to be considerably larger than n% Note that in our approach 
also the onset of direct J/ip suppression will be abrupt, governed by the cluster size at 
ric{ ss (J/ip). - Since the very loosely bound ip' can in addition also be dissociated by 
hadronic comovers we use below n c cl the experimental suppression for the 8 % of the 
J/i/j from this state. 

We now calculate the distribution of nucleon-nucleon collisions and the corresponding 
cluster density and cluster size in S — U and Pb — Pb collisions as function of the impact 
parameter b. In the Table, we show the results for r = 0.20, 0.26 and 0.30 fm. It is seen 
that the cluster density for central (6 = 0) S — U collisions is reached in Pb — Pb collisions 
for an impact parameter 8.0fm<6<8.5 fm. This result is found to be quite independent 
of r, even though the associated cluster densities and sizes show a strong r-dependence. 
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14.1 
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7.62 
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149.0 


878.4 


605.2 


179.9 


136.5 


40.2 


A cl 
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0.30 


23.7 


111.0 


80.0 


28.2 


22.0 
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n d /a 


(fm- 2 ) 




6.32 


7.91 


7.57 


6.40 


6.22 


5.60 



Table: Collision densities and cluster sizes for S — U and Pb — Pb collisions. 



Tuning r in the given range leads to variations in the size of the cluster at the onset 
of deconfinement, and this in turn determines how abruptly the J/ip production rate 
drops at that point. We find best agreement with the most recent data [^5| for r = 0.26 
fm. With a total nuclear overlap area of 35 - 40 fm 2 for both central S — U collisions 
and Pb — Pb collisions at b = 8 fm, the cluster size at the onset of deconfinement is for 
this r value less than half the total area; hence the density profile varies little within a 



given cluster. The resulting form for the anomalous J/ip suppression is shown in Fig. 3. 
As direct consequence of the finite cluster size at the onset of deconfinement we get an 
abrupt onset of the suppression as function of impact parameter b. As function of the 
measured transverse energy Et, the drop is softened by the inherent Et — b smearing 



of the experiment, which is here included in the form described in [18]. The agreement 
between our model and the data is seen to be quite good. 
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Figure 3: The anomalous J/tp suppression given by our model, compared to data |25|j ; 



the righthand figure shows the results after division by pre-resonance absorption |18 



Finally we return to the relation of deconfinement and percolation. For an infinite 
two-dimensional continuum system, percolation occurs at n = 1.175. For r = 0.25 fm 2 , 



we obtain n c ~ 1.2, for r = 0.26 fm 



1.3; thus the threshold value obtained from 



central S — U collisions is quite close to the infinite area percolation point, supporting the 
conjecture of percolation as basis for deconfinement [@-@, [11 . 



In conclusion: the onset of deconfinement for a sufficient overlap of strings appears to 
provide a conceptually clear picture for the beginning of quark-gluon plasma formation. 
Its main and model-independent consequence is that a critical collision density implies 
clusters of a critical size. The resulting abrupt onset of deconfinement is found to agree 
with recent data on anomalous J / ip suppression in Pb — Pb collisions. 
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